        program OISM1
C Ocean model small mode number 1: Reff=0.1mic sigma=0.4	
        character*10 ipa1
        real xtau,n(20),k(20),wldisc(20)
        real rad(22),dvdlnr(22)
        real cij1,cij2,rm1,rm2,rst1,rst2
	real pi,cte,tmp1,tmp2
      data wldisc /0.350,0.400,0.412,0.443,0.470,0.488,0.515,0.550,
     s             0.590,0.633,0.670,0.694,0.760,0.860,1.240,1.536,
     s             1.650,1.950,2.250,3.750/
        data n/1.45,1.45,1.45,1.45,1.45,1.45,1.45,1.45,
     s         1.45,1.45,1.45,1.45,
     s         1.45,1.45,1.43,1.43,1.43,1.43,1.43,1.43/	
        data k/0.0035,0.0035,0.0035,0.0035,0.0035,
     s         0.0035,0.0035,0.0035,0.0035,0.0035,	
     s         0.0035,0.0035,0.0035,0.0035,0.0035,	
     s         0.0035,0.0035,0.0035,0.0035,0.0035/	
        data rad /0.050000,0.065604,0.086077,0.112939,0.148184,
     &  0.194429,0.255105,0.334716,0.439173,0.576227,0.756052,
     &  0.991996,1.301571,1.707757,2.240702,2.939966,3.857452,
     &  5.061260,6.640745,8.713145,11.432290,15.000000/
        pi=atan(1.)*4.0
	cte=sqrt(2*pi)
	rm1=1.48
	rst1=0.60
	cij1=1.0
	write(6,*) "22"
	do i=1,22
	tmp1=exp(-(log(rad(i))-log(rm1))**2./2./rst1/rst1)
	dvdlnr(i)=cij1*tmp1/cte/rst1
	write(6,*) rad(i),dvdlnr(i)
	enddo
	write(6,101) (n(i),i=1,20)
	write(6,102) (k(i),i=1,20)
 101    format(20(F10.5,1X))
 102    format(20(F10.7,1X))	
	stop
	end
	 
       

